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"Solar  Fterct  and  MafurtoaplMric  Partidca: 
iBvcstigatioiis  Based  Upon  the  ONR-d(tt  and  ONR-604  Experiments" 


I.  Introduction: 

This  performance  rqiort  covers  work  acoxi^lished  under  ONR  Grant  N00014-90-J-1466 
related  to  the  radiatkxioivitonment  in  near-Earth  space.  The  goal  of  tiie  imeardi  is  to  measure  and 
describe,  quantitatively,  die  Geospace  radiati(»  environment  in  which  men  and  ^lacecraft  must 
survive  and  function.  The  ttxdsfcMr  this  investigation  are  die  data  returned  by  die  ONR-602  and 
ONR-604  experiments,  both  flown  uiuler  the  au^ices  of  ONR  and  the  Air  Force  Space  Test 
Program,  augmented  by  correlative  databases  of  bodi  ^[lace-based  and  ground-based  data.  The 
investigation  involves  data  analysis,  modeling  and  (^plications  to  a  variety  of  space  equipment  and 
environments. 

This  repot  builds  upon  the  detailed  Tedinical  Report  (Fall,  1992)  and  the  previous 
perfomance  repcxts.  For  the  current  period,  the  principal  effort  was  in  the  modeling  of  the 
Anomalous  Con^xxient,  determining  the  modulatiai  level  during  die  CRRES  Misskm,  and 
studying  die  quiet-time  particle  access  to  the  magnetosphere. 


II.  The  Anomalous  Conqxinent: 


As  described  in  the  previous  report,  it  is  necessary  to  have  a  representation  of  the  Anomalous 
Component  (AQ,  particularly  the  anomalous  helium  in  order  to  undmtand  the  solar  modulation 
level  as  a  function  of  time.  Iherefore,  we  have  undertaken  die  task  to  develop  realistic  local 
interstellar  spectra  (LIS)  for  the  anomalous  cosmic  rays. 

An  anomalous  con^xment  to  the  observed  cosmic  ray  spectrum  has  been  reported  for  the 
elements  helium,  carbon,  nitrogen,  oxygen,  neon,  and  argtxi.  All  six  of  these  elerronts  have  high 
first  ionization  potentials  and  are  presoit  in  the  interstellar  medium  as  neutral  atomic  matter. 
Amxnalous  component  cosmic  rays  are  widely  believed  to  originate  as  neutral  atoms  in  interstdlar 
space,  which  drift  into  the  inner  solar  system  and  are  subsequently  photoionized  by  solar  UV 
radiation  and  swept  baidt  outward  by  the  solar  win  .  Since  the  AC  species  are  ionized  by  solar 
UV,  one  expects  Ae  AC  cosmic  rays  to  be  singly  charged. 


Acceleration  at  the  shock  boundary  of  the  heliosphere  produces  a  local  "interstellar" 
spectrum,  which  is  then  modulated  in  the  heUosphere  as  die  energetic  particles  move  back  toward 
Earth.  The  energies  ccxisistent  with  acceleration  at  the  helio^heric  shock  boundary  are  low,  and 
the  AC  local  interstellar  spectra  are  consequendy  very  steep,,  essentially  vanishing  at  energies  above 
about  1  GeV  per  nuclecm.  The  low  energy  of  dw  AC  cosmic  rays  accounts  fin-  tteir  near  absence 
from  measund  cosmic  ray  fluxes  during  periods  of  high  solar  modulation.  However,  the  AC 
contribution  at  low  energies  has  been  observed  in  both  the  mid-70's  and  the  mid-80's,  during 
successive  miniina  in  the  solar  cycle,  and  is  expected  to  return  during  each  solar  minimum  poiod. 


Fcht  our  pre^t  woric,  it  was  decided  to  develqi  AC  local  interstellar  spectra  (xily  for  helium, 
nitrogen,  oxygra,  and  neon.  The  observed  AC  fluxes  for  carbon  and  especially  i^<m  are  quite 
low,  and  the  inclusion  of  anomalous  components  oi  these  species  would  not  significant  Tte  AC 
ccmtribution  is  more  si^uficant  for  the  remaining  four  elements,  and  in  particular,  it  is  important  to 
include  the  AC  helium  in  the  analysis  of  dw  solar  modulation  level  as  a  functicm  time,  as 

described  in  the  noct  section.  For  each  of  the  four  elements  included,  die  modd  AC  spectrum 
contains  only  the  primary  isotqie  (Le.,  ^He,  l^N,  1®0,  and  20Ne).  Contributions  other 
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isotc^)es  would  be  veiy  small  cooqpared  to  the  primary  isotope,  and  ctxisidering  that  die  AC  flux 
itself  is  conqiaiatively  small  exc^  at  the  lowest  ener^es,  the  added  refinement  of  attempting  to 
model  the  isou^ic  composition  is  not  warranted  at  this  stage. 

The  helium  and  oxygen  local  interstellar  spectra  were  derived  at  die  University  of  Chicago, 
by  demodulation  of  measmed  data  (Garcia-Munoz,  private  ctxnnaunication)  and  have  been  ad^ted 
for  this  analysis.  The  LIS  for  nitrogen  and  neon  were  obtained  by  renormalizing  the  oxygen 
spectrum,  with  the  normalization  factcnr  chosen  to  yield  die  best  fit  obtainaUe  to  a  variety  of 
measured  AC  datasets  corresponding  to  different  tevels  of  modulation  and  different  locations  in  die 
helio^here.  It  was  possible  to  obtain  reasonable  fits  to  the  datasets  for  these  diree  elements  at 
energies  above  about  10  MeVAiucletxi  using  just  a  single  spectral  shape,  that  of  the  demodulated 
oxygen  spectrum.  Measurements  fircxn  IMP-8  at  1 AU,  the  Voyager  spacecraft  at  2  and  at  21  AU, 
and  Pioneer  10  at  16  AU  were  used  in  fitting  the  AC  spectra.  Details  of  these  fits  are  summarized 
in  Table  1. 

Figure  1  shows  the  Local  Interstellar  Spectra  (US)  that  have  been  adopted  for  the  four 
anomalous  elements  that  have  been  studied.  The  LIS  are  cut-off  at  S  MeVAiucleon,  but  the  N,  O 
and  Ne  would  be  expeaed  to  turn-over,  like  the  He,  at  lower  energies.  Modeling  this  turnover  has 
not  been  possible  within  the  available  data.  Figures  2-9  show  the  results  of  our  predictions,  using 
the  spectra  of  Figure  1  as  input,  compared  to  the  measurements  cited  in  Table  1.  The  LIS  were 
mod^ated  to  the  time  interval  and  the  location  in  the  heliosphere  of  the  spacecraft  fic»n  which  the 
data  were  reported.  The  agreement  between  the  prolictions  (solid  lines)  and  the  data  points  is 
reasonably  good  except  for  Voyager  H  results  from  1985-86  at  21  AU  fw  N  and  O  below  10 
MeVAiucleon.  The  equivalent  data  for  anomalous  Ne,  however,  is  in  good  agreement  with  the 
predictitxis.  Whether  this  indicates  a  problem  witii  the  N  and  O  data  during  1985-86  or  with  the 
predictions  is  unknown.  The  agreement  for  E>10  MeV/nucleon  is,  however,  good  enough  for  the 
purposes  of  the  present  study. 

Note  that  the  AC  is  singly  charged  and  has,  therefore,  a  higher  rigidity  than  a  Galactic 
Cosmic  Ray  (GCR)  ion  of  the  same  energy.  In  order  to  combine  the  AC  spectra  with  the  GCR 
spectra,  a  routine  was  written  to  generate  a  binary  file  that  mimics  the  output  of  the  GCR 
prc^gation  calculation.  The  solar  modulation  jxogram  accqits  this  file,  modulates  die  elonents 
with  a  charge  of +1,  then  outputs  the  spectra  on  the  same  energy  grid  as  the  GCR  spectra.  The 
modulated  AC  spectra  can  tlwn  be  added  point-by-point  to  the  G^  spectra,  if  desi^.  It  is 
important  to  note  that  the  AC  must  be  modulated  s^arately  fiom  the  GCTl  component,  because  of 
its  singly-charged  icmization  state. 

During  the  CRRES  period,  the  highest  level  of  modulation  occurred  at  the  beginning  of  the 
missitm.  Taking  the  period  fipcxn  8/19  - 12/15/90,  the  mean  modulation  level  is  characterized  by 
value  of  <|>  =  1589  MV.  Figure  10  shows  spectra  for  He,  N,  O,  and  Ne  modulated  to  this  level. 
The  solid  curve  includes  the  AC  while  the  dashed  curve  does  not  As  expected  from  previous 
measurements,  there  is  essentially  no  AC  helium  present  A  surprise,  however,  is  the  N,  O  and 
Ne  which  each  show  a  small  amount  of  AC  just  above  10  MeV/nucleon  even  at  these  high  levels  of 
modulation.  The  explanation,  we  believe,  is  that  being  singly-charged,  the  rigidity  for  N,  O  and 
Ne  is  so  much  larger  than  for  He  that  the  modulation  cannot  stqp  some  penetration  of  these  itxis  to 
1  AU  even  at  solar  maximum!  Understanding  this  effect  quantitatively  \rill  require  additiraial 
analysis  (unfortunately  ONR-604  cannot  clwck  the  result  because  it  is  not  sensitive  at  these  low 
energies),  but  this  may  also  help  to  explain  tli^  source  of  the  "trapped  anomalous  oxygen"  repented 
recently  from  measurements  on  the  SAMPEX  spacecraft  (Chimmings  et  al.,  GRL,  20, 2003, 

1993). 
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Spectrum: 
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Hefinm 

Demodulated  at  UnivemycfCJiicago 

See  reference 

Garda-Munoe.  et  aL 
21st  ICRC,^L  164,1990 


Oxygen 

Demodulated  at  UniveniQr  of  Chicago 

IMP-8. 1  AU,  1977,  Phi »  S06  MV 
Voyager  I  fe  n,  2  AU,  9/77  -  2/78,  Phi »  S13  MV 
Honcer  10. 15  AU.  l/H  -  3/78,  Phi «  513  MV 
Voyager  n.  21  AU,  10/85  -  8/86,  Phi »  693  MV 

Gaicia-Munoz,  et  aL,  unpuUished  figure 

Cummings,  Stone,  and  Webber,  Ap  J,  287.  L99  -  L103, 15  Dec.  1984 

Cuinmings  and  Stone,  Solar  Wind  VI  Proc.,  NCAR/TN  -  306f  PROC 
V.  J.  Pizzo,  Ed..  1988 


Nitrogen 

Normalized  point-by-pmnt  fiom  anomalous  oxygen  qtectrum: 

Nitrogen  flux  =  (Oxygen  flux)/4.502 

Voyager  I  &  n,  2  AU,  9/n  -  2/78,  Phi  s  513  MV 
Voyager  n.  21  AU.  10/85  -  8/86.  Phi  =  693  MV 

Cummings,  Stone,  and  Webber,  Ap  J,  2SL  L99  -  L103, 15  Dec.  1984 

Cummings  and  Stone,  Solar  Wind  VI  Proc.,  NCAR/TN  -  306+  PROC, 
V.  J.  Hzzo,  Ed.,  1988 


Necm 

Normalized  pmnt-by-pdnt  fiom  anomalous  oxygen  plectrum: 

Neon  flux  » (Oxygoi  flnx)/47.147 

Voyager  I A  n.  2  AU,  9/77  -  2/78,  Phi  =  513  MV 
Voyager  11,21  AU,  10/85  -  8/86,  Phi »  693  MV 

Cummings,  Stone,  and  Webber,  Ap  J,  287.  L99  -  L103, 15  Dec.  1984 

Cummings  and  Stone,  Solar  Wind  VI  Proc.,  NCAR/TN  -  306+  PROC, 
V.  J.  Pizzo,  Ed.,  1988 


Flux  (m'^  sr  s  MeV/Nucleon) 


Flux  (m*^  sr  s  MeV/Nucleon) 


Fig.  2.  Anomalous  Oxygen  at  1  AU 
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Fi£.  4.  Anomalous  Oxygen  at  15  AU 


Fig.  5.  Anomalous  Oxygen  at  21  AU 
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Fig.  7.  Anomalous  Nitrogen  at  21  AU 
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III.  Solar  Modulation  Level: 


In  the  previous  report  we  described  an  evaluation  of  the  level  of  solar  modulaticm  ~  the  0 
parameter  —  based  upon  the  IMP-8  helium  flux.  That  evaluation  suffered  from  the  presence  of 
anomalous  helium  in  the  data,  which  was  not  included  in  the  calculations.  Employing  the  final 
GCR  helium  spectra  and  the  anomalous  helium  component,  described  above,  we  have  re-done  the 
analysis.  The  result  is  shown  in  Figure  1 1,  extending  frran  late  1973  through  early  1993,  the  full 
extent  of  the  available  data.  The  points  represent  measurements  of  ^  matte  by  other  groups  and 
extracted  from  the  literature  or  from  private  conomunicadons.  Now,  the  agreement  is  excellent! 

We  believe  that  this  comparison  fully  validates  the  solar  modulation  analysis  and  lends  confidence 
to  the  use  of  the  derived  ^  parameter  in  the  interpretation  of  the  data. 

One  of  our  objectives  from  the  study  of  the  radiation  environment  in  Geospace  has  been  to 
use  the  available  data  and  the  models  to  predict  the  environment  that  will  be  encountered  by 
spacecraft  We  can  proceed  toward  that  goal  for  the  GCR  plus  AC  if  we  can  predict  the 
irKxlulation  parameter  ^  (MV).  In  order  to  predict  future  levels  of  modulation,  Ae  assumption  was 
made  that  the  modulation  level  will  be  periodic,  following  the  22-year  solar  cycle  as  discussed  in 
the  previous  report.  The  values  of  <|>  for  the  interval  January  1970  through  May  1971  are, 
therefore,  assumed  to  reflect  the  values  from  January  1992  through  May  1993.  Likewise,  (|)  from 
November  1995  through  December  2010  can  be  obtained  frrom  the  corresponding  values  for 
November  1973  (when  the  IMP-8  helium  data  began)  through  December  1988.  Most  of  the  40 
year  interval  1970-2010  can  thus  be  filled  in  using  the  helium-based  phi-versus-time  profile  in 
Figure  11. 

In  order  to  avoid  predicting  exactly  the  same  transient  features  seen  in  the  phi-versus-time 
plot,  the  set  of  monthly  average  <!>  values  from  1973-93  was  smoothed  using  a  Fast  Fourier 
Transform  method,  and  the  smoothed  set  of  values  was  then  translated  22  years  into  the  past  and 
22  years  into  the  future  to  provide  estimated  monthly  average  values.  Figure  12  (top)  shows  these 
intervals,  together  with  the  established  values  from  1973-93.  As  is  seen  in  the  figure,  there  remain 
two  gaps  of  identical  length,  from  June  1971  through  October  1973  and  from  June  1^3  through 
OctobCT  1995,  reflecting  the  fact  that  the  IMP-8  data  available  at  present  fall  29  months  shon  of 
spaiming  a  full  22-year  solar  cycle. 

As  we  have  shown  previously,  there  is  a  good  correlation  between  the  modulation  level  and 
the  Climax  (or  any  other)  neutron  monitor  counting  rate,  for  which  data  exist  back  to  the  1950's. 
Therefore,  we  us^  the  neutron  monitor  for  the  interval  6/71-10/73  as  a  guide  in  filling  the  gap. 

We  assume  that  the  overall  qualitative  behavior  of  the  neutron  monitor  and  the  modulation  will  be 
similar  over  any  time  interval  not  chosen  too  short  Even  over  spans  of  monAs,  there  are  changes 
in  neutron  monitor  rate  that  are  not  immediately  mirrored  in  the  helium  fluxes,  and  vice  versa. 
Therefore,  the  neutron  monitor  infoimation  was  used  only  as  a  check  to  help  ensure  that  the  <|> 
values  estimated  to  fill  the  gaps  do  not  fail  to  reflect  any  major  trends.  The  neutron  monitor  data 
was  snooothed  using  the  same  method  mentioned  previously,  so  that  just  the  qualitative  trend  of  the 
neutron  count  rate  can  be  seen.  Monthly  values  of  the  modulation  level  for  6^1-10/73  and  6/93- 
10/95  were  then  fitted  into  the  gaps  by  hand,  attempting  more  or  less  to  reproduce  the  trend  of  the 
neutron  monitor  data  and  avoid  any  discontinuities  in  the  modulation  level  at  either  end  of  the  gaps. 

The  resulting  total  set  of  monthly  average  levels  of  nxxlulation  for  the  entire  interval  1970- 
2010,  with  both  remaining  gaps  filled,  is  shown  in  Figure  12  (bottom).  This  is  now  the  database 
used  for  calculating  average  values  of  the  level  of  nxxlulation  over  any  specified  time  interval.  The 
accuracy  of  these  predictions  will  be  checked  only  by  future  spacecraft  measurements. 
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Hgurell.  The  derived  moduladcm  parameter  ^(MV)  as  a  functioned  time  comptttd  to 
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Hgure  12.  Moclhilatkm  parameter  ^vemis  time  shoiwing  the  ptosected  values  based  iq)m  a  22  ynrsdar 
cycle  (top)  and  the  full  result,  1970  -  2010,  with  the  gaps  filled  as  described  in  the  text 
(bottom). 


rv.  Particle  Pen^tkm  into  Magnetoq)hae: 


Part  of  our  effcxt  over  the  past  quarter  has  been  devoted  to  omtinuing  the  study  of  particle 
access  into  the  inner  regions  of  the  magneto^here.  We  have  been  focusing  on  Helium  since  this  is 
a  high  abundance  elonent  and  may  be  relatively  straight-ftxward  to  understand.  The  raw  data, 
excluding  sdar  active  periods,  is  shown  in  Figure  13  as  the  live-time  normalized  helium  flux 
versus  L  shell.  Direct  access  for  Galactic  Cosmic  Rays  is  observed  from  high-L  values  witfi  a 
rapid  decrease  in  the  region  of  L  =  3-4  representing  the  cut-off  of  the  GCR  access.  What  is  most 
interesting  are  the  peaks  at  L  ~2  and  L  ~1.3  well  inside  dte  region  availaUe  to  GCR  events.  These 
low-L  helium  are  probably  a  "tnq)ped"  component. 

The  helium  events  are  derived  from  the  priority  P2  pulse-height  analyzed  events  after  "clean¬ 
up",  i.e.  consistency  checks,  cuts  on  the  penetration  flags  and  angle  cuts,  as  described  in  a 
previous  report  The  resulting  tracks  are  shown  for  L>6  in  Hgure  14  where  the  AE-E  matrix  is 
plotted  for  both  ID=7  (stq)ping  in  Kl,  preceding  detectt»’  is  D6, 1  mm  thick)  and  ID>7  (stoi^ing 
in  K2-K8,  preceding  detector  is  5  mm  ^ck).  In  both  cases,  there  is  excellent  charge  resolution 
with  the  helium  well  separated  from  the  Li  and  heavier  ions.  Q4ote  diat  protons  are  analyzed  in 
priority  P3,  so  are  not  included  in  this  P2  analysis  except  for  a  few  background  events  which  ate 
eliminated  by  the  "clean-up"  procedure.) 

Using  the  trajectory  system  that  is  included  in  the  ONR-604  instrument  on  CRRES,  we  can 
measure  the  arrival  direction  of  each  particle.  Knowing  the  spacecraft  location,  attitude  aird  spin 
period  we  can  projea  the  arrival  direcdcxi  onto  the  local  magnetic  field  directitm  to  determine  the 
particle's  local  pitch  angle,  a.  The  result,  shown  in  Figure  15  for  the  helium  with  L  <  2.3,  is  a 
distribution  sharply  peaked  near  cos  a = 0,  i.e.  at  a  values  near  90*  or  270*.  Such  a  distribution  is 
characteristic  of  a  trapped  particle  pq)ulation  and  demonstrates  that  the  low-L  events  are,  at  least, 
pseudo-trapped  in  the  magnetosphm. 

This  leads  directly  to  the  question  of  the  time  dependence  of  the  intensity  of  this  trapped 
component  The  CRR^  results  are  shown  in  Figure  16  (top)  compared  to  the  Climax  neutron 
monitor  counting  rate  (bottom).  The  Climax  data  give  an  idea  of  the  state  of  the  helio^)here  by 
recording  the  high  energy  cosmic  ray  intensity.  During  the  first  part  of  the  (31RES  mission,  the 
neutrcm  nwnitor  rate  increased  slowly  as  die  solar  cycle  moved  from  solar  maximum  towards  solar 
minimum.  During  this  period  the  trapped  helium  intensity  decreased  by  ~30%.  This  argues  that 
the  low-L  helium  (a)  is  not  being  replenished  by  galactic  cosmic  rays,  and  (b)  is  a  slowly  decaying 
population. 

The  cosmic  ray  recovery  was  interrupted  by  the  major  solar  events  of  March-June  1991 
ir<cluding  the  3/24/91  substtxm  which  fonned  a  "new"  radiation  belt,  as  reported  from  CARRES 
data.  During  this  period  the  trapped  helium  continued  to  decrease.  This  shows  that  the  helium  was 
not  injected  by  the  1991  SEP  events  as  had  been  suggested  previously.  If  anything,  it  could  be 
argued  (withii  d.e  very  small  statistics)  that  some  of  this  tr^ied  hdium  was  reii»ved  by  the 
substorm.  Throughout  the  rest  of  the  CRRES  period  the  trapped  helium  continues  to  be  present  at 
a  reduced  intoisity. 

The  suggesticMi  of  a  solar  flare  injection  origin  cannot  be  totally  mled  out  by  Figure  16. 
Rather  we  can  say  that  the  1991  flares  did  not  jnovide  any  injection.  Whether  a  flare  such  as  tiiose 
at  the  end  of  1989  could  have  provided  the  particles  (and  we  are  observing  the  decay  of  this 
component)  cannot  be  determined  from  die  present  dataset 


We  have  divided  die  data  into  two  grcaips  at  odMt  S85,  the  peak  cf  die  flare  March,  1991. 
The  normalized  L  dependence  for  these  two  periods  is  shown  in  the  bottom  two  panels  ci 
Figure  17.  Note  that  before  the  March  flare  events  at  L~1.9  are  dcxninant  while  after  the  flare 
(even  considering  the  small  stadsdcs)  the  L~1.9  peak  is  almost  gone.  Thus,  the  decrease  in  die 
trapped  helium  rate  shown  in  Figure  16  appears  to  be  due  almost  entirely  to  a  reduction  in  die  peak 
at  L^l.9.  A  more  detailed  time  history  be  required  to  investigate  diis  effect  in  detail 

The  distrilMitions  in  Hgure  17  may  be  die  edges  of  a  single  distribution,  centered  upon  the 
inner  radiation  belt  For  diis  to  be  the  case,  the  ONR-604  instrument  would  have  had  to  be 
insensitive  to  events  occurring  between  L~1.3  and  L~1.8.  This  could  result  if  the  pulse-hei^t 
analysis  was  disabled  due  to  high  detector  counting  rates,  so  that  the  logic  circuits  never  received 
the  appropriate  signals  for  P2  analysis.  To  begin  to  look  at  this  possibility,  the  top  panel  of  Hgure 
17  plots  the  P3  counting  rate  as  a  Unction  of  L.  (Note  die  logarithmic  scale.)  Above  "lO^  in  this 
rate,  we  could  expect  problems  in  the  logic  and  die  inner  edges  of  the  Helium  distributions 
correspond  to  rates  of  a  few  times  10^.  Thus,  it  is  possible  that  what  we  are  observing  is  in  fact  a 
region  of  trapped  high  energy  helium  co-located  with  the  inner  proton  radiation  belt  but  slightly 
broader  in  L  dian  the  proton  belt. 

Qearly  additional  analysis  will  be  required  to  verify  these  tentative  conclusions.  This  will  be 
one  of  the  thrusts  of  our  work  for  the  next  quarter. 


CRRES  Live  Time  Normalized  Helium  Flux  at  Various  L  Regions 


Hgure  13.  The  normalized  L  distribution  of  quiet-time  Helium  evmits. 
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I^ue  14.  P2  dita  in  AE>E  format  for  He<>fioinqiiiet-tiines.  Results  are  shown  for  1D»7  (top) 
and]D>7  (bottom)  for  orbital  inMcvalswidiL>6. 
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